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INTRODUCTION 

This Study describes the simulation on the EAI 
HYDAC® 2400 HYbrid Digital-Analog Computer 
of a control system for a tubular chemical reac- 
tor. Specifically, it describes the transient simu- 
lation of the reactor itself, its product separator, 
and its control system (Figure 1 shows a block 
diagram of the physical system simulated); the 
study calls for the solution of distributed param- 
eter reactor equations, which can be accomplished 
best by a hybrid computer having the mathematical 
capability of function storage. 

Although the chemical process industries (CPI) 
have made extensive use of the analog computer 
alone in the past to solve complex design, control, 
and data analysis problems, the simulation of dis- 
tributed parameter systems impose aheavyburden 
on the computing capability of the hardware usually 
available. In most instances, for example, the 
capacity of a normal computer facility (say 100 - 200 
amplifiers) rarely is challenged by problems re- 
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Figure 1: Block Diagram of Physical System 



quiring the simulation of linear and non-linear 
ordinary differential equations. 

However, mathematical models involving several 
partial differential equations (PDE) ...models of 
fixed bed transients, tubular reactors, etc fre- 
quently tax this computing capacity, and the avail- 
able hardware definitely becomes the limiting 
factor in such cases, usually requiring simplifica- 
tion of the equations, the overall system model, 
or both, and some compromise of the simulation 
accuracy. These excessive hardware requirements 
for the solution of PDE' s are due to the transforma- 
tion of each PDE into a set of ordinary differential 
equations (ODE), each of which must be solved 
simultaneously (i.e., in parallel). This requires the 
implementation by means of analog components of 
a large number of almost identical computing cir- 
cuits, a requirement which is usually sufficient 
to exhaust the equipment complement of even the 
most fully supplied laboratory. 

This situation has been changed by using a different 
mathematical formulation of the equations which 
requires a function memory capability. Hybrid 
computers, coupling analog speed with digital 
memory and accuracy have satisfied this require- 
ment and thus enable system analysts to utilize 
the maximum problem solving capability of their 
computing facilities. In addition, they have pro- 
vided simple means for implementing logical op- 
erations, transport delays, and the transient 
simulation of stage-wise processes... processes 
which are, themselves, hybrid in nature. 

The purpose of this Study is to illustrate the solu- 
tion of a typical chemical process problem on the 
HYDAC 2400 Hybrid Computing System, consisting 
of an EAI 231R-V General Purpose Analog Com- 
puter, a DOS-350 Digital Operations System, and a 
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DCS-375 Digital Computing System. The problem 
described is the simulation of a tubular reactor 
(represented by 2 PDE's), its separator, and a 
control system. Problem definition, organization, 
programming, and hardware requirements are pre- 
sented in detail. 

OBJECTIVES 

The objectives of this simulation are typical of 
those encountered in many CPI simulations. They 
are: 



c , using a control system to vary the feed rate of 
component C. Since B is present in excess, vari- 
ations in C (note reaction IB) have a direct effect 
on product formation. The control system selected 
is a conventional three-mode controller. 

It has been assumed that the separator is 100% 
efficient and that its dynamics can be represented 
by a first-order lag. Sensor and valve dynamics 
associated with the control loop are assumed 
negligible and the valve-position-versus-flow re- 
lationship is considered to be linear. 



(1) To define the physical system mathematically, 

(2) To propose a control system to maintain a 
constant output product flow rate, and 

(3) To determine controller settings for the pro- 
posed control system. 

PHYSICAL SYSTEM* 

The complete system under study, shown schemat- 
ically in Figure 1, consists of a tubular reactor, a 
separator, and a control loop. The kinetics of the 
system are: 



A + B 



C + D 



(1A) 



MATHEMATICAL MODEL 

Assumptions: The assumptions used to establish 
the mathematical model of the problem are as 
follows: 

1. The velocity through the reactor, V, is 
constant. 

2. The process is isothermal with the con- 
sequence that the rate of reaction is 
dependent only on composition. 

3. The density of the reaction mixture, p, 
and the average molecular weight of the 
reaction mixture, m, are constant, inde- 
pendent of composition. 



C + B 



E + D 



(IB) 



where component E, the desired product, is formed 
by feeding A, B, and C to the reactor. Component 
D is a waste product which is discarded. 

The reactor effluent is fed to a separator, which 
ideally separates the effluent into five pure product 
streams, Aj Ej. Unreacted A and C are re- 
cycled and combined with make-up feeds, Ai and 
Ci, to obtain reactor feed rates, A Q and C Q . Com- 
ponent B is returned to tankage on leaving the 
separator. It is supplied in excess to the reactor 
at a constant feed-rate, B G . Excessive amounts of 
B are used to maintain favorable reactor behavior. 

Disturbances (step changes and noise such as valve 
chatter or low frequency sinusoids) enter the sys- 
tem through the make-up feed stream, Aj, and 
cause variations in product flow, Ex. The purpose 
of the control loop is to operate on the flow error, 



4. The reactor is of the plug-flow type; i.e., 
there is no axial diffusion. 

The first assumption can be justified, even though 
the flow rates of A, B, and C at the inlet to the 
reactor vary, if one of the following conditions 
exists: 

1. The flow of one component, say B, is very 
high (and constant) compared with the 
combined flow rate of the other compon- 
ents. This is the case for the problem 
under consideration. 

2. All the reacting species (A, B, C, D, E) 
are dissolved in an inert carrier which 
forms the bulk of the flowing stream. 

Kinetics: The rates of conversion for the compon- 
ents are given by: 



r = -k„ ab + k.cd 
a 1 3 



(2a) 



r fe = .-k ab + kgCd - k 2 bc + k 4 de (2b) 



'This system was patterned after the one described in Reference (1). 



r = k, ab - k„cd - kJbc + k.de 
c 1 3 2 4 



(2c) 



2 



r. =k,ab - k n cd + kJx> - k.de 
d 1 3 2 4 

r = k rt bc - k.de 
e 2 4 



(2d) 
(2e) 



where a e = concentrations of components A, 

...E, lb moles/ft 
r ,...,r = rates of formation, lb moles/ 

a 6 (ft 3 ) (sec) 
k l'"*' k 4 = reaction rate constants, (ft 3 )/ 
(lb mole) (sec) 

Only two of these equations are required, since 
r b, r c> anc * r d can be stoichiometricly related to 
r and r . 



where 

A E = the flow rates of each component, 

moles/ sec 
7 = cross- sectional (or flow) area, ft 
F = volumetric flow rate, ft 3 / sec 
p = fluid density, lbs/ft 3 
m = fluid molecular weight, lbs/lb mole 

Since A = Fa and B = Fb, equations 4 and 5 can 
be rewritten in terms of flow rates to obtain 



v — = k — +k — 

at Sx IF 3 F 



(7a) 



r, = r - r 
b a e 



r = -r - r 

c a e 



r = -r 
d r b 



Integrating these equations yields 



b=b + a - a - e 
o o 



c=c + a -a-e 
o o 



d = b - b 



(3a) 
(3b) 
(3c) 

(4a) 
(4b) 
(4c) 



The initial concentration, d , is zero and is omitted 

o 

from equation 4c. 



il*? + V — = k — DE 
Bt 3x 2 F " 4 F 



B = B - A + A - E 

o o 

C=C +A -A-E 

o o 

D = B - B 
o 



(7b) 

(7c) 
(7d) 
(7e) 



Separator Equations: The separator equations, 
which are first-order lags, are, in transfer func- 
tion notation. 



E. 



(8) 



A L °L 



E tS+ 1 
L 



Reactor Equations: The concentration dependence of 
components A and E with respect to time and posi- 
tion is expressed by two partial differential 
equations: 



where 



-1 



S = Laplace operator, sec 
t = separator time constant, sec 



da , TT Sa 
— + v = r 

at a x' a 



ae a e 

— + v — = r 

at ax e 



(5a) 
(5b) 



where 



t = time, sec 

x = distance from entrance of reactor, ft. 
V = velocity of mixture through the reactor, 
ft/sec 

Assuming V, P, andm are constant, z (total molar) 
flow rate also is constant; therefore, 



The subscripts "L" and "1" denote reactor and 
separator effluent concentrations, respectively. 

Separator equations were not written for compon- 
ents B and D since their effluent concentrations 
are not required in the simulation. 

Control System: The classical transfer function for 
a 3-mode controller is 




r S + l 
r 



aT D S + 1 /Vr S + K 



(9) 



where 



z = £S = A + B + C + D + E 



m 



Fp 

— = CONSTANT (6) 



= derivative time constant, sec 

T = reset time constant, see 
r 

7 = derivative gain, dim 



3 



= reset gain, dim 

K = controller gain, sec/lb mole 
c 



A = A, + A + f 
o 1 s 



/C \A1 

Al = change in percent of valve opening, dim _ „ n [ m ) 

o = C l + s + \100/ 



(15a) 



(15b) 



The controller gain includes the gain of the valve, 
sensor, valve positioner, etc. The constants K c , 
r r , and are settings available to plant opera- 
tors whose determination is one of the objectives 
of this study. 



The error signal, e, is 
e = E - E, 



(10) 



where E s , the set point, is the desired product flow 
rate. 

Feed and Recycle Equations: Based on the linear 
valve assumption (as shown in Figure 2) , the make- 
up feed, Cj, of "C" is 



C. = C + 
l s 



Al 



(11) 



where 



C = design flow rate of make-up feed, moles/ 
sec 

C = maximum flow rate of make-up feed, 
m 

moles/sec 



The make-up feed rate, A of component "A" is 

A. = A + f (12) 
is 

where 



A = design flow rate of make-up feed, moles/ 
sec 

f = disturbance introduced into the system, 
moles/sec 

The disturbance can take two forms: a step change, 
f = K (13) 

or a sinusord, 

f = K Sin cot (14) 

where 

K = magnitude of the disturbance, moles/ sec 
60 = frequency of the disturbance, rad/sec 

Equations 11 thru 14 can be combined with recycle 
flows, A^ and C to obtain reactor feed equations. 



Recall the feed rate, B , of component "B" is 
fixed. 



Method of Characteristics: The partial differential 
equations can be reduced to ordinary differential 
equations by the method of characteristics (2). Con- 
sider the first order partial differential equation 



a F SF 

If + V (x, t) - = g (F. x, t) 



(16) 



which is of the same form as equations 7a and 7b. 
The solution is wanted for t>0 and 0<x<L, where 
L represents the length of the reactor. 

Throughout this general discussion, F will be used 
to describe the flows (A,...,E) in the mathematical 
model. 

The method of characteristics requires that inte- 
gration be performed along lines in the space-time 
domain, defined by 



(17) 



To find the equivalent of equation 16, write the total 
derivative of F 



at 3x 
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to obtain 



or 



dF dF 3F dx SF SF _ xX % 
dt = Ft + B7dt = si +V ^(F,x,t) (18a) 



dF a f a F 
V d^ = 3T + V 37=S(F,K,t) (18b) 



This equation is an ordinary differential equation 
and can be solved directly by an analog computer. 
For the problem under consideration, V is constant 
which means that the characteristics consist of 
straight lines having a slope of V, as shown in 
Figure 3. The initial condition, F (0, t), is given by 
the problem, and for the tubular reactor, this 
represents the variation in the flow of components 
at the inlet to the reactor. Since the conditions at 
the outlet of the reactor are desired, the procedure 
based on the method of characteristics is to inte- 
grate equation 19b on the computer until x = L, or 




o L 
Figure 3: Characteristic Lines in the Space-Time Domain 

For example, if F (0, t 0 ) represents the condition 

at the reactor inlet at time t Q , the integrator which 

generates F would be reset to F (0, t G ), and then 

allowed to integrate until x = L, at which time the 

integrator would produce the values of F at the 

outlet of the reactor t seconds later, i.e., F (L, 

t + t). 
o 

Since the characteristic curves for the problem 
under investigation are straight lines having a slope 
of V, the residence time in the reactor is a con- 
stant (which also is the computation time, L/V). 



High speed computation techniques will be used to 
allow a large number of characteristic lines in 
the space-time domain to be "scanned" in order 
to obtain a semi-continuous solution. 

Computational Model: The reactor- separator sys- 
tem is presumed to be operating at its design point 
when a disturbance is introduced into the system. 
Therefore, good scaling practice dictates that equa- 
tions be rewritten in perturbation form. Defining 



AA L 


= A L 


A Ls 


(19a) 


AC L 


= C L 


" C Ls 


(19b) 


AE L 


= E L 


Ls 


(19c) 


AA 1 


= A 1 


- A ls 


(19d) 


AC 1 


= C 1 


- C ls 


(19e) 


AE 1 


= E 1 


- E = -e 

s 


(19f) 



yields the following separator-feed equations 
AC. 



A A, 



AC, 



A E 



TS + 1 



A = A + A A, + f 
os 1 



C = C 
o s 



+ AC 1 + ( 



C \A1 
m i 

100/ 



(20a) 
(20b) 
(20c) 



where 



s and C s denote A + A - ~ 

S lS S -LS 



and C + C. 



respectively. 



The non-linearities in the reactor equations make 
perturbation variables impractical; however, their 
form lends itself to solution by the method of 
characteristics. Applying this technique (see equa- 
tion 18b) yields 



^= kl AB-k 3 CD 



= j^- = k DE - k BC 
dy 4 2 



(21a) 
(21b) 



where 



y = 



_ y x 



X 

vf 



The above equations, solved in conjunction with the 
three algebraic equations (equations 7c thru 7e) 
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defining B, C, and D in terms of A and E, are the 
computational model of the reactor. 

The controller equation (equation 9) requires no 
modification since it is already in perturbation 
form. However, it should be noted that Al is 

limited (o<l = 1 + A 1 <100). 

s 

As mentioned before, a high-speed computational 
procedure is required to allow the reactor outlet 
product flows to be obtained in a semi- continuous 
manner. 

METHOD OF SIMULATION 

High Speed Computing Techniques: Referring to the 

description of the Method of Characteristics and to 
Figure 4, if integration along the characteristic 
line is speeded up so that it occurs in a small inter- 
val of time, A (including reset), the value of F at 
the end of this interval would be F (L, t 0 + A). If 

this value is delayed by — -A seconds and then in- 
troduced into the simulation, real-time simulation 
of F (L, t + ~ ) is accomplished. In the HYDAC 

2400 system, the obvious means for storage is to 
use the digital computer PCS 375) memory. 




-a -4 — V-a— ! 



Figure 4: Diagram of High-Speed Computing Technique 

The benefit of high-speed computation results 
from the fact that over the time inerval L/V, a 
large number of computations (namely L/VA) can 
be made instead of just one computation being 
performed in real-time. This number can be as 
high as several hundred in a practical case. It has 
to be chosen equal to the number of memory words 
available in the delay simulator(3) divided by the 
number of data points to be stored. In the present 



case, 150 memory words were used for each flow 
rate variable of interest. 

Computational Assignments: Table 1 indicates the 

computational assignment of each part of the 
HYDAC 2400 System. Note that the analog com- 
puter operates in two speeds. It simulates the 
separator, controller, and feed equations in real- 
time while using high-speed repetitive operation 
to solve the reactor equations. 



Table 1: Computation Assignments in the 
HYDAC 2400 System 



Unit 


Function 




High-speed solution of Equa. (21) 
for simulating tubular reactor. 


231-R 
(Analog) 


Real-time simulation of separator 
and controller. 


Problem control including: 

1. Control of repetitive operation 

2. Control of DOS modes 


DOS-350 


Control of information between DOS 
and System 375 through use of digital- 
analog and analog -digital converters . 


SYSTEM 

375 
(Digital) 


Storage of data for the delay required 
to simulate the tubular reactor. 



The computer circuits to be described are based 
on the requirement that three of the flow rates of 
components (A, C, and E) are to be simulated in 
real-time. This can be seen from Figure 1 in 
which the reactor outlet flows of only these com- 
ponents affect the inlet flows (i.e., A andC are re- 
cycled and E is used to adjust the supply of fresh 
C). This means that the delay simulator in the 
digital part of the computer must be able to handle 
three channels of information. 

The retention time in the reactor, (L/V), is to be 
3.0 sec. If the operate period is 12 ms and the 
reset period is 8 ms, then the combined period, A, 
is 20 ms. The number of storage locations in the 
DCS 375 computer for each variable (A, C, and E) 
is 3000/20 or 150. 

Problem Control Assignment: Since an investigation 
of this kind requires continual modification of both 
forcing functions (f) and parameters (controller 
settings), and continuous function readout, etc., 
problem control has been assigned to the analog 
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computer. Its mode of operation controls both the 
DOS -3 50 System and DCS 375 computer. 

COMPUTER PROGRAMMING 

Analog Computer Programming: The scaling, pro- 
gramming, and problem check-out of the analog 
portion of this study is straightforward and will not 
be discussed at this time. However, a table of 
scaled voltage equations is shown in Table 2 for 
information purposes, and Appendix A contains cir- 
cuit diagrams, pot sheets, etc. 

Referring to the computer diagram, track- store 
units 12, 13, 17, 18, 22, and 23 monitor the flows 
of A, C, and E and stored AA L *, ACl*, and AE L *, 
the terminal flow rates. This occurred during the 
operate period (which was 12 ms) of the high- speed 
repetitive operation (HSRO) cycle when y=l. It was 
implemented using the logic signal from an elec- 
tronic comparator (ECOO). The terminal values 
were held constant during the 8 ms reset period 
of the HSRO cycle and transferred to digital 
memory. The integrators in the HSRO portion of 
the simulation were speeded up by a factor of 100. 

Table 2: Summary of Scaled Voltage Equations 



REACTOR 

d [lOOy] 
~d^ 

-d [50A] 
dr 

-d [50E] 
dr 



(»[f -] • OH 

[25D]= 2.5^- I ^[l00] - [lOB]j 

[10B] = (^)[100] -(i)[ 5 0 (A o - A)] -(i)[ 5 0E] 
[50C] = |50 cJJ + |~50 <A q - A)j - [50E] 



SEPARATOR 



[ 50 AA J "> f 50AC ]] 10 
[5AA L ]"- S + i: [5AC L ]" TS + 1 

[50 c] -10 

[ 5 ae l ]"- +1 



CONTROL LOOP t 

[2Al] [l0 €l ] [2c 2 ] [2Al] 
- x x 

[lOOeJ [lOOcJ [10 ^e 2 ] 
[lOeJ S + 1/t d 

[iooc] 

[ 2g 2] 
[10 c J 

[2AlJ 



(10 0!) S + 10/ T. 



D 




S + 1/t t 



[2€ 2 ] 
FEED EQUATIONS 



S + 1/K t 
r r 



[soa o ] = [so aaJ + (^)[ioo] + [sof] 

o] ■ t AC,] +(%00] ^f 4 j 



|50C 



DISTURBANCE EQUATIONS 
Step Input 



Sinusoidal Input 



[50f] - (f)[l00] 
put 

[50f] = Sin to ,t] 



Notes: ( ) = pot setting, [ ] = computer variable or reference volt- 
age, and T— machine time = j8y in reactor equations* 

t Since L is assumed to be 50%, 2 A ? is limited to ± 700 volts. 



Track-store units 27 and 28 held A Q and C Q con- 
stant during the operate portion of a HSRO run. 
During this period their input was updated, but the 
updated values were not used in the present com- 
putation. They represented new feed flows for the 
next "slug" of fluid integrated down the reactor 
length. 

All digital to analog conversion (DAC) was per- 
formed during the operate portion of the HSRO 
cycle, while analog to digital conversion (ADC) 
was performed during the reset period. 

Integrators involved in the "real-time" simulation 
of the separator and control system have a unity 
R-C time constant. 

Signals generated on the MLG patch panel used in 
the DOS program include the mode of computer 
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operation (OP, HD, IC). A logical one appeared 
at the mode termination when the computer was in 
a particular mode. The HSRO, A, and A" signals 
were required for conversion purposes. The A 
signal was "high" during the reset period and 
"low" during the operate period (A is the comple- 
ment of A) . 

The pulse train which produces A and A operates 
continuously regardless of the computer mode; 
therefore, it must be synchronized with the DOS 
clock if it is to be used on the DOS to initiate 
conversion operations. 

DOS Programming: The programming of the DOS 
is best discussed by considering the principle 
areas of the over-all program, which is shown in 
Appendix A. These areas are ADC, DAC, and con- 
version timing; the initialization abstracts of the 
over-all circuit diagram will be used in the dis- 
cussion of the over-all DOS program. 



ADC 



The circuit shown in Figure 5 is designed to con- 
vert three analog voltages, appearing on trunks 
100-102, after receiving a blip input signal. A blip 
completion signal is generated after the three con- 
versions are completed. 

The blip input resets the multiplexer, selects chan- 
nel zero, and simultaneously calls for a conver- 
sion (CONV input to ADC). When the ADC is com- 
plete, a FIN signal is gated with a ready signal 
from the DCS 375 Computer (OC 12) to enable 
(TE^ the transfer of information to the digital 
memory. When this occurs, the Input-Output Buffer 
Register (IOBR) between the DOS- 3 50 and the 
DCS 375 is simultaneously cleared and loaded in a 



parallel manner (Lp, CLR inputs to IOBR). The 
ADC is reset (RST) as is OC 12 at this time. 

When transfer into memory is complete a blip ap- 
pears at TCj which increments (INC) the multi- 
plexer to channel one and initiates the conversion 
cycle. After the third conversion has been com- 
pleted a blip appears at the output of AND gate 28B 
signifying completion of all three conversions. 



DAC 



Consider Figure 6, where a blip input sets a flip- 
flop (211) to signify conversion is required and 
simultaneously clears (CLR) the general purpose 
counter (GC 10) which is a channel selector for the 
DAC. A digital ready signal (OC 11) through an 
"AND" gate (240) enables transfer out of the digi- 
tal computer into the IOBR. A blip from the TC G 
termination loads the digital word into the assem- 
bly register (L p ), whose completion is signified 
by a LFIN-RST blip which initiates conversion 
(TRA) to channel zero. When conversion is com- 
plete a blip from T FIN- RST advances the BCD 
counter to channel one as well as initiating another 
DAC conversion by sending a reset signal to flip- 
flop 220. 

After three conversions have been completed, the 
AND of the TFIN-RST signal with channel 2 signal 
(AND gate 28A) produces a blip finish signal, which 
resets flip-flop 211 inhibiting further conversion 
signals. 



Con 



version Timing 



The OP, HD, IC, and A signals received from the 
MLG unit (as shown in Figure 7) are passed 
through AND gates as a safeguard against noise 



IOBR 


Lp 


S 


1 ■ 




■12 


CLR 




231 Rj DOS 



RST 



s 


1 ■ 




-12 


FIN 


ADC 




0 


1 


2 


3 


CONV 



RST 
f0CI2| 



I RST 



8 


4 


2 


1 


MXPR 



INC 



OUT 
-O 



RST 



-O 
IN 



Figure 5: ADC Circuit 
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or spikes. The OP and A signal are sent to flip- 
flops to insure that they are synchronized with the 
DOS clock. 

Temporarily disregarding initialization, after an 
ADC is complete (flip-flops 201 is high) and the 
rep. op. cycle is in operate (flip-flop 221 is low), 
a blip is sent to the DAC program through the 
initializing circuit. After the DAC program is 
complete, flip-flop 210 is set and waits for the 
reset part of the rep. op. cycle. At this time a blip 
appears at the output of AND gate 21A to initiate 
the ADC routine. The ADC completion signal sets 
flip-flop 201 and the above procedure is repeated. 



This is accomplished by DIF 004 (Figure 7), which 
is passed through the OR circuit feeding the DAC. 
After the first set of DAC is complete and the start 
of the next reset period is reached, an ADC is 
initiated. DIF 003 is used to prevent an ADC from 
taking place after the reset period has started in 
the event that the first DAC was completed during 
a reset period. 

The digital program, which will be discussed later, 
replaces values removed by DAC with the next 
ADC values and then advances to the next set of 
memory locations. 

Digital Computer Programming 



Initialization 

To insure proper operation of the system, the cor- 
rect initial values of reactor recycles must appear 
prior to the first ADC. If the reactor feed voltages 
are not correct, its effluent voltages will be in- 
correct. In setting up the problem it has been as- 
sumed that prior to an entering disturbance the 
reactor is operating at steady-state. Therefore, 
the digital memory cells used for the delay are 
filled with the steady- state effluent flow rate 
values. In lieu of the above, a DAC must take place 
prior to an ADC. 



The stored-program digital computer program 
provides storage and playback of the AA^, ACj^, 
and AE L reactor outlet flow rates. The total delay 
required is 3 seconds and each high-speed analog 
run requires 20 ms.; therefore, a total of 150 
points must be stored per variable. Since the re- 
actor is assumed to start from steady- state, the 
digital program also loads the table initially with 
the steady-state values (in this case zero). 

The flow diagram for the digital program is shown 
in Figure 8 and the assembler program is listed 
in Appendix A. 



IA 



^NIT^- 



ENABLE 
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Figure 8: Flow Diagram of Digital Program 
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The computer is started from location INIT where 
interrupts are enabled and the interrupt trap cells 
are loaded with the appropriate jump instructions. 
The program then goes to STRT where the initial 
steady- state values are loaded into the table, and 
the initial address of the table is established in 
memory. Table indexing is then initialized after 
which the program proceeds to the D to A trans- 
fer routine. The index is then stored in memory 
for use in the A to D transfer routine. 

The first data word then is loaded into the ac- 
cumulator and the parallel output channel is en- 
abled (OC 11 is set to indicate to the DOS that the 
DCS 375 is ready to output a word). The output 
channel is tested for the ready flip-flop set, which 
halts the digital program until the DOS commands 
the data to be transferred. When the channel be- 
comes ready, the DCS 375 outputs the word to the 
IOBR of the appropriate DAC channel. 

The index, N, is tested then for readout of the last 
DAC channel. If it is the last channel, the program 
proceeds to prepare for the A to D conversions. 
If not, N and i are indexed and a test is performed; 
i is tested for the end of the table. At the end, an 
error halt will occur since the routine should have 
exited at the N test. If the index is non-zero, the 
routine loads the accumulator with the next word 
and repeats the output cycle. 

The A to D transfer program first establishes the 
proper index, i, from the memory location QF and 
reinitializes the variable count, N. The parallel 
input channel then is enabled and OC 12 is set to 
inform the DOS that the program is ready to input 
analog data. The program then enters awaiting loop 
until the input channel ready signal is brought high 
at the completion of an A to D conversion. 

The program inputs the data word from the IOBR 
and stores it in the proper indexed location. The 
variable count, N, is tested for 3 and, if it is not 
the last channel, the routine increments N and i 
and tests i for the last table address. At the last 
table address, the program executes an error halt 
since an exit should have occurred at the N test. 
For a non-zero, the program returns to point 5 
to prepare to input the next channel. 

For a variable count N =3, the routine reinitializes 
N, increments i, and tests for the last data point. 
If it is not the last point, the program returns to 
point 2 to prepare for the next D to A transfer. If 
the index is zero, i is reinitialized before re- 
entering the D to A transfer routine. 

Computer Control: When the analog computer is in 
the IC mode the DOS is in CLR; it can only be put 



into the RUN mode if the analog is in operate, OP, 
or hold, HD. If the DOS is put in RUN, an interrupt 
signal is fed to the digital computer, thru I^> 
which causes the digital program to clear its 
memory and load delay cells with initial, steady- 
state values of A, C, and E. DIF 002 disables the 
interrupt signal until the analog computer again 
goes from the IC to the HD mode. 

This computer control set-up uses the analog com- 
puter as the master computer. It permits param- 
eters to be varied on the analog without clearing 
the DOS or reprogramming the digital computer. 
To obtain this flexibility of operation, the HD mode 
cannot be by-passed; to start a run, one must go 
through the HD mode. 

Since both the ADC and DAC are completed in less 
than 1 ms, there is no need to remain in HD for any 
predetermined period of time. 

COMPUTER RESULTS 

The results of the simulation are shown in Figures 
9 thru 15 and are summarized in Table 3. 

Runs 1 through 3 (Figures 9 thru 11) are illustra- 
tive. They show typical steady- state reactor con- 
centration profiles, the L/V transport delay in the 
reactor, and the insensitivity of the control vari- 
able to a 10-radian-per-second sinusoid. This 
proves that the separator will filter out high fre- 
quency disturbances. 

Run 2 proves that the step input disturbance does 
not result in instability, but establishes a new set 
of steady- state operating conditions. The function 
of the control system, of course, is to maintain the 
design outflow of E in spite of this change. 

Adding proportional control (run 4), indicates that 
a gain of 10 will yield a steady- state off- set whose 
magnitude is much less than that obtained without 
proportional control. 

Increasing the gain to 20 (run 5), results in peri- 
odic response (period equal to 9 sec) indicating 
that the system is at or approaching its critical 
gain (-180° phase shift). Applying the Ziegler- 
Nichols^ criteria for a proportional plus integral 
controller resulted in approximate controller set- 
tings of K =10 and t =5 seconds. 



tK c = 0.45 (CRITICAL GAIN), ? r = 0.8 (PERIOD OF OSCILLATION 
AT CRITICAL GAIN) - in this case, K c = 70, r.= 5 sec 
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Figure 9: 



Steady-State Concentrations Versus Reactor 
Length Variable 



System response to these settings, show in Figure 
7, is reasonable, and the controlled variable re- 
turns to its set-point in about 20 sec. Further re- 
finements could have been attempted but these 
results were adequate and satisfied the objectives 
of the study. 

It should be noted that the response obtained in run 
7 is typical of that expected for a load change im- 
posed on a controlled system. 

Run 6 justifies undertaking a study of this type 
since it proves that unstable control settings exist. 
Failure to maintain proper control can of course 
lead to: 



equations were required in the simulation, a con- 
ventional all-analog approach would have required 
two complete 231-R computers. Only one analog 
computer would be required in the hybrid simulation 
of the system if the thermal equations were added. 
Therefore, one advantage of the hybrid computer 
simulation of partial differential equations is a 
definite saving of analog components over the con- 
ventional parallel approach. 

Table 3: Summary of Computer Results 
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1 


Typical steady- state concentration 
versus reactor length 


10 


2 


System response to a step change 
(no control system) 


11 


3 


System response to a sinusoidal 
disturbance (no control system) 


12 


4 


System response to a step change 
with proportional control; K c =10 


13 


5 


System response to a step change 
with proportional control, K c =20 


14 


6 


Example of system instability with 
improper control settings 


15 


7 


System response to a step change 
with proportional plus integral con- 
trol; K = 10, r =5 
c ' r 



1. Increased operating costs due to high re- 
cycle rates, low production of product, etc. 

2. Potential danger to operating personnel 
due to toxic materials, fires, etc., which 
could result from a run- away or im- 
properly controlled system, 

which are typical motivations for applying computa- 
tional techniques to chemical process problems. 

The advantages of simulating this system on a 
hybrid computer are worthy of comment. A parallel 
simulation of the reactor equations would have re- 
quired a complete 120-amplifier, 231-R computer 
to perform this simulation; the analog requirement 
for the hybrid simulation was less than 50 ampli- 
fiers. If the reactions were exothermic and thermal 
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Figure 10: Run 1, System Response to Step Input 
(No Control System) 




Figure 11: Run 3, System Response to Sinusoidal Input 
(No Control System) 




Figure 12: Run 4, System Response to Step Input 
(K c = 10) 




Figure 13: Run 5, System Response to Step Input 
(K c = 20) 




Figure 14: Run 6, System Behavior with Improper 
Control Settings 




Figure 15: Run 7, System Response to Step Input 
(K c = 10, Tr = 5) 
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